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Abstract. We extract constraints on the pion distribution amplitude from available data on the pion- 
photon transition form factor in the framework of light-cone sum rules. A pronounced discrepancy 
(2.7 — 3) C between the Gegenbauer expansion coefficients extracted from the CELLO, CLEO, 
and Belle experimental data relative to those from BaBar is found. Predictions for the pion-photon 
transition form factor are presented by employing a pion distribution amplitude obtained long ago 
from QCD sum rules with nonlocal condensates. These predictions comply with the Belle data but 
Q_i disagree with those of BaBar beyond 9 GeV 2 . 
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PION-PHOTON TRANSITION FORM FACTOR IN LCSRS 

The pion-photon transition y* + y* — > n° (see [1] for a review) is defined by the correlator 
of two electromagnetic currents 

J d 4 ze-^(n°(P) | r{j M (z) Jv (0)} | 0) = ie^ vap qfq P 2 F^((f, q 2 ) , (D 

where F^* 7 * 7t (Q 2 ,q 2 ) is the pion-photon transition form factor (TFF) with the photon 
momenta q\ and q 2 , and Q 2 = — q\ > 0, q 2 = —q\ > 0. This transition process was 
measured by the CELLO [2], CLEO [3], BaBar [4], and Belle [5] Collaborations for 
the kinematics Q 2 ^> mp, q 2 <C nv 2 ^ (the same also for the planed measurement by 
BESIII [6]). For this kinematics, perturbative QCD factorization is reliable only in the 
leading-twist approximation, while higher twists are also important as shown in [7]. For 
q 2 < m 2 , one has to include the interaction of the (quasi) real photon at long distances 

~ 0(1/a/?). To this end, we apply [8, 9] 2 Light Cone Sum Rules (LCSR)s [18, 19] 
that effectively account for long-distance effects of the real photon making use of quark- 



1 Presented by the first author at the QCD@Work Conference, Lecce (Italy), 18-21 June, 2012. 

2 This contribution is based on our previous works [8, 9], where the detailed description of our LCSR 
application and the comparison of our results with some of other pion-photon TFF studies [10-17] can be 
found. Here, we only briefly describe the basic ingredients of our LCSR approach. 



hadron duality in the vector channel and applying a dispersion relation in q 2 to get [19] 
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where so ~ 1.5 GeV 2 is the effective threshold in the vector channel and M 2 is the 
Borel parameter (0.7 — 0.9) GeV 2 . The effect related to the small nonzero momentum 
q 2 was recently discussed in the framework of LCSRs [20], and also using Monte Carlo 
simulations [21]. Here, we neglect this effect and apply the real-photon limit q 2 — > that 
can easily be taken in the LCSRs given by (2): Fyf^Q 2 ,q 2 ) = i r yy*^(2 2 ,0). 

The spectral density has been calculated in QCD up to the level of the twist-six (Tw-6) 
term [15] (using similar notations for the other twists): 
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where the various twist contributions are given in the form of a convolution of the hard 
parts with the pion distribution amplitude (DA) of a given twist. At the same time, for 
the twist-two contribution we have the perturbative expansion 
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where a s = a s /4iz. In our studies all calculated terms have been included and they are 
shown in Fig. 1 as a chain of ovals each enclosing an area proportional to the absolute 
value of the corresponding contribution. The sum of all terms of each row gives its total 
contribution, as detailed in the accompanying table in this figure. 
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FIGURE 1. Relative percentage contributions to the pion-photon TFF Fyy< K (Q 2 ) with the asymptotic 
pion DA at Q 2 = 3, 8, and 25 GeV 2 in table form and graphically (shown only for Q 2 — 8 GeV 2 ). The 
blue color corresponds to negative terms, while the red one denotes the positive terms. 

As we see from Fig. 1, the twist- six contribution [15] and the next-to-next-to-leading 
order (NNLO) term [22] proportional to the pVpart have absolute values of similar 



size but opposite signs in the whole <2 2 -region. Therefore, we consider their sum as 
a contribution to the theoretical uncertainties both for fitting the pion DA and for 
estimating the pion TFF for a given pion DA model. All other terms are included into 
the central value of the TFF: Leading order (LO), next-to-leading order (NLO) of twist- 
two [22] (corrected in [15]), and twist-four [19]. The uncertainties of the applied LCSR 
approach contain the sum of the NNLO twist-two and twist-six contributions, a 20% 
variation of the twist-four coupling 5 2 (/i 2 ) = 0.19 GeV 2 [23], and those uncertainties 
stemming from the pion DA model. 

We consider the pion DA in terms of the first two coefficients 3 of the full Gegenbauer 
expansion (Jc= 1 — x): 

(p K (x,n 2 ) = 6xx 1 + a2(l~L 2 )C^ 2 (x - x) + a4(n 2 )C^ 2 (x - x) + ... 

The coefficients a n (jl 2 ) are evolved according to the Efremov-Radyushkin-Brodsky- 
Lepage (ERBL) [24, 25] evolution equation to the NLO level [26]. 



Data analysis and predictions for the TFF 

The described LCSR-based techniques allow the computation of the TFF for any pion 
DA given in the form of a Gegenbauer expansion. Using LCSRs with pion DAs obtained 
in QCD sum rules with nonlocal condensates [27], one obtains the predictions shown by 
a shaded (green) band in Fig. 2 in comparison with the available experimental data. 
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FIGURE 2. Theoretical predictions extracted in [8, 9] for the scaled y* yx TFF, calculated in the LCSR 
approach with pion DAs obtained in QCD sum rules with nonlocal condensates [27], in comparison with 
the experimental data from CELLO [2], CLEO [3], BaBar [4], and Belle [5]. 

On the other hand, one can fit the pion DA Gegenbauer coefficients to the different 
data sets in order to obtain corresponding experimental constraints. This can also reveal 



The consideration of the pion in terms of three coefficients can be found in [8, 9]. 



the extent of the discrepancy among the various data sets. In Fig. 3, we present the con- 
fidential region of and given in the form of error ellipses derived from fitting dif- 
ferent sets of data. Best-fit values of the ^-squared goodness of fit j} /ndf (ndf=number 
of degrees of freedom) are shown at the centers of the ellipses, while deviations between 
particular data sets are displayed on the sides of the triangle in units of one standard 
deviation (la ~ 68%). From Fig. 3, we conclude that the inclusion of the BaBar data 
to those of CELLO&CLEO leads to a 3 a shifting of the confidential region from the 
(blue) ellipse at the bottom on the right to the (red) solid-lined ellipse at the top on the 
left, accompanied by a significant increase of the ^-squared goodness of fit value from 
0.4 to 2. The additional inclusion of the Belle data plays no role in the two-parameter 
analyses (dashed-dotted (red) ellipse). If we include to the CELLO&CLEO data only 
the Belle data (ignoring the BaBar data), then the shifting of the confidential region is 
moderate (1 .2a), with a slight increase of the ^-squared goodness of fit value from 0.4 to 
0.6 (dotted black ellipse). Note that we considered the ^-squared goodness of fit for the 
analyzed experiments, ^ 2 , as the sum of the individual ^-squared goodness of fit values 
associated with each experiment. When the data strongly deviate from each other, this 
way of combining data could be inefficient for extracting the confidential region of the 
fitting parameters, though it is still good for identifying and rating this deviation. 





^4-rr T - :r ^_+BaBar 






0.1 
















0. 




2.7o- 


col 






BELLE | 


0.1 








' 3.0<r\ — U - - - r>* 








1 1.2a ~ 




0.2 




# 0.4 








\. CELLO & CLEO 




0.3 














^Cl2 



0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 



FIGURE 3. Confidential regions for two-parameter {a%, 04) fits based on three data sets: first - CLEO, 
CELLO [2, 3] (solid (blue) ellipse at the right bottom); second - CLEO, CELLO, and BaBar [4] (solid 
(red) ellipse at the left top); third - CLEO, CELLO, and Belle [5] (dotted ellipse at the center). Best fit 
values of the ^-squared goodness of fit x 1 jndf are shown at the apexes of a triangle. On its sides we 
display the discrepancy among the data sets in terms of a standard deviation (la fts 68%). The area of the 
a%, 04 values allowed by QCD sum rules with nonlocal condensates from [27] is indicated by a slanted 
(green) rectangle. 

From Fig. 3, we conclude that the BaBar data deviate from the other data sets at 
the 3a level and cannot be satisfied by pion DAs based on models with only two 
Gegenbauer coefficients — unlike all other data sets — and hence they require (at least) 
a sizable coefficient or still higher coefficients. 



CONCLUSIONS 



It is not a trivial matter that QCD has come to use collinear factorization to describe basic 
exclusive processes, like the pion-photon transition form factor. This scheme has been 
challenged by the BaBar data because they indicate an auxetic behavior at Q 2 > 9 GeV 2 
that cannot be uniquely described by a known QCD mechanism. Therefore, the Belle 
data, which do not replicate this increase of the scaled pion-photon transition form factor, 
may help restoring the confidence to the standard QCD scheme, though a definitive 
answer requires more accurate data in the high-<2 2 range. 
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